Validity of the rigid band picture for the t—J model 

R. Eder, Y. Ohta and T. Shimozato 
Department of Applied Physics, Nagoya University, Nagoya 464-01, Japan 

We present an exact diagonalization study of the doping dependence of the single 
particle Green's function in 16, 18 and 20-site clusters of t—J model. We find evi- 
dence for the validity of the rigid-band picture starting from the half-filled case: upon 
doping, the topmost states of the quasiparticle band observed in the photoemisson 
spectrum at half-filling cross the chemical potential and reappear as the lowermost 
states of the inverse photoemission spectrum. Features in the inverse photoemission 
spectra, which are inconsistent with the rigid band picture, are shown to originate 
from the nontrivial point group symmetry of the ground state with two holes, which 
enforces different selection rules than at half-filling. Deviations from rigid band be- 
haviour which lead to the formation of a 'large Fermi surface' in the momentum 
distribution are found to occur at energies far from the chemical potential. A Lut- 
tinger Fermi surface and a nearest neighbor hopping band do not exist. 
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A well known problem in the description of high-temperature superconductors is the 
volume of the Fermi surface. Since these systems are close to a metal-to-insulator transition, 
there arises the following question: should one model them by a system of quasiparticles 
which correspond to the doped holes and populate the dispersion relation calculated for a 
single hole (rigid band approximation) or should one assume that the ground state can still 
be obtained by adiabatic continuation from the noninteracting one, so that the Fermi surface 
corresponds to a slightly less than half-filled band of noninteracting electrons? Based on 
numerical studies of the momentum distribution and single particle spectral function for the 
frequently used t—J model it has been argued ^ that the single hole represents a 'problem 
of only marginal relevance' for the doped case: already for two holes in clusters with 16—20 
lattice sites (corresponding to a nominal hole concentration of ~10 %) a kind of phase 
transition has ocurred so that both, the Fermi surface and the quasiparticle band structure in 
its neighborhood, resemble that for noninteracting particles. In this manuscript we present 
evidence against this widely accepted picture: as far as the photoemission spectrum is 
concerned the rigid band approximation (RBA) is in fact an excellent one; to fully understand 
the inverse photoemission spectrum one has to take into account the nontrivial symmetry 
of the two-hole ground state, which enforces transitions into a second, symmetry-different 
band of many-body states. 
The t—J model reads 

H = -t + H.c.) + J E (^^ ■ - 

<i,j>,a <i,j> 

The Si are the electronic spin operators and the sum over < i,j > stands for a summation 
over all pairs of nearest neighbors. The operators q^o- are expressed in terms of ordinary 
fermion operators as Ci^„{l—ni_^). We study the single particle spectral function An{k, u) = 
An-{k, —u) + A„ _|_(fc, cu), where the photoemission (PES) spectrum An-{k,Lj) and inverse 
photoemission (IPES) spectrum An,+ {k,uj) are defined as 
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Here \'^u,n) (E^n) is the z/*'' eigenstate (eigenenergy) with n holes (in particular z/=0 implies 
the ground state) and the reference energy Eref is chosen as i?o,n (although other values will 
be advantageous later on). All spectra were evaluated exactly by the Lanczos method. 
Let us first recall a few constraints due to sum rules. The momentum distribution 
n(fe) = (\E'o,n|Cfc o.Cfc,o-|^o,n) in the n-hole ground state is the zero*^ moment of An-{k,u!); 
the zero*'' moment of An^+{k,uj), is given by l—N/2L—n{k), where N is the number of 
electrons and L the number of sites (for an even number of electrons and total z-spin 0). In- 
troducing Q={n, tt) and 6n{k)=n{k)—n{k + Q), the expectation value of the kinetic energy 
can be written as 



where e{k) is the free-particle energy, and the summation over k is restricted to the interior 
of the antiferromagnetic Brillouin zone (AFBZ). Since e(fc)<0 in the range of summation, 
in order to have {Ht)<0 the average of 6n{k) in the AFBZ must be positive i.e. n{k) must 
be larger inside the AFBZ than outside, so that there is always a tendency towards the 
formation of a 'large Fermi surface' in n{k). To exemplify this. Tab. |, compares n{k) for 
the two-hole ground state of the fully isotropic t—J model and of the t—Jz model (where 
the transverse part of the Heisenberg exchange is discarded) with an added staggered field. 
The latter model has explicitely broken symmetry (the ground state expectation value of the 
staggered magnetization is 61.3% of the value for the Neel state) which rigorously excludes 
a large Fermi surface. In spite of this, n{k) is almost indistinguishable for the two models 
and in particular always would suggest a large Fermi surface. This shows that use of n{k) to 
assign a Fermi surface may be quite problematic. The key problem is, that for strongly 
correlated systems the 'quasiparticle peak' near the chemical potential carries only a small 
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fraction of the total PES weight; if we decompose n{k) into a component from the coherent 
peak, nQp{k), and from the integration over the incoherent continua, riindk), we usually 
have ninc{k) > nQp{k). Then, when only n{k) is considered, for example a change in 
ninc{k) (which is unrelated to low-energy physics) may mimick a Fermi surface. As will be 
seen below, this is precisely what happens when the t—J model is doped with holes. Next, 
the sum rules for the spectral function enforce that in the doped case PES weight (IPES 
weight) is concentrated inside (ouside) of the AFBZ, so that not only n{k) but also the 
distribution of spectral weight in fc-space is reminiscent of free electrons. 
Let us next discuss what can be reasonably expected if the RBA were to hold: a variety of 
diagonalization studies pj-Q have shown that two holes in the t—J model form a bound state 
with a binding energy (Eb^O-SJ—J) that is a sizeable fraction of the single-hole bandwidth 
(iy~2J). The two-hole ground state thus should be modelled by a state of the type |$o) = 
X]fc A(fc)a][ ^Ofc Jo), where a^,^ is the creation operator for a hole-like 'quasiparticle' in the 
band observed in AQ _{k). The wave function A(fc) may differ appreciably from zero for all 
quasiparticle states within r^Ep above the ground state so that a 'Fermi surface' does not 
exist. In IPES we annihilate a hole (with momentum k) so the remaining hole should be 
in a state belonging to the single hole band with momentum —k; the intensity should be 
proportional to the quasiparticle occupation n(fc). Conversely, in photoemission we add a 
third hole and, neglecting all interactions except for the Pauli principle, should observe the 
same spectrum as for a single hole with the intensity of the peaks near the Fermi energy Ef 
being reduced by a factor l—n{k). Adding up the weights of the peaks in A2^+{k,uj) and 
A2-{k,u!) closest to Ef we should therefore ideally obtain the weight of the 'unsplit' peak 
in y4o,-.(fc, cu), a sensitive quantitative test for the RBA. Moreover, the dx2_y2 symmetry of 
the two-hole ground state 0] implies that in the 16 and 18 site cluster A(A;) must have a 
node along the (1, 1) direction and hence fi(fc)=0 for these momenta. 

Let us now check, in how much these expectations are borne out by the exact spectra. Fig. 
|T] compares A(fc, uj) in the half-filled and two-hole ground state for all allowed momenta in 
the 4x4 and 18-site cluster ((0, 0) and (vr, vr) are from the 18-site cluster). Fig. || displays the 
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same information for the 20-site cluster. No adjustment of any kind has been performed. 
In agreement with the above discussion Fig. |1| shows that for momenta along the (1, 1) 
direction (left panel) there is a striking similarity between the PES spectra for the doped 
and undoped case, A2^-{k,uj) and y4o.-(fc, cj), particularly near the Fermi energy i?/~+1.7t. 
A band of peaks with practically identical dispersion and weight can be clearly identified 
in both groups of spectra (a possible exception is this momentum might play a 

special role since the single hole ground states at this momentum and (vr, vr) have total spin 
S = 3/2, whereas all other single hole ground states are dubletts [^). Away from (1,1) 
(right panel of Fig. |l]) doping leads to a shift of weight from the PES band to IPES peaks 
immediately above Ef-. this suggests the 'split peak' situation n{k)j^O. The situation is the 
same for the 20-site cluster. Fig. |^ shows the same similarity between the low energy parts of 
the PES spectra for the doped and undoped case. Again the the dominant low energy PES 
peaks near the Fermi energy (right panel) remain either unaffected by doping or seem to 
(partially) cross the chemical potential to reappear as low energy IPES peaks. Point group 
symmetry poses no significant constraint on the hypothetical pair wave function A (A;) in 
this cluster, (all momenta except (0, 0) and (vr, vr) have low symmetry) and consequently we 
observe low energy IPES weight for all momenta where low energy PES peaks were present 
at half-filling. 

It is only at energies remote from Ef that major changes in the PES spectra occur upon 
doping: for momenta around (0,0) and (tt, vr) there is a pronounced addition/depletion of 
incoherent spectral weight at energies ~3t below Ef. The increase/decrease of the integrated 
spectral weight near (0, 0) / (vr, vr) and the corresponding formation of a 'large Fermi surface' 
in n{k) upon doping therefore is clearly not the consequence of the phase transition-like 
emerging of a nearest neighbor hopping band [|I| in the range 2 J ^ 3t around Ef (note that 
the spectral weight near Ef at (0, 0)/(7r, it) even diminishes/increases upon doping). Rather 
it is accomplished by the reshuffling of incoherent spectral weight deep below Ef, and thus 
is certainly unrelated to any low-energy physics. As a matter of fact, the actual form of 
A{k,Lj) for the doped 20-site cluster rules out the Luttinger Fermi surface postulated ^ 
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for this cluster on the basis of the momentum distribution. Despite the fact that (vr, 0) 
and (7r/5,37r/5) are on opposite sides of the Luttinger Fermi surface, the low energy peak 
structure in A{k,uj) is practically identical for these momenta: both photoemisson spectra 
show a pronounced low energy peak, indicating that these momenta are 'occupied'. The 
criterion n{k) < 1/2 (n(fc) > 1/2), employed in Ref. to decide whether a fc-point is inside 
(outside) the Fermi surface, obviously has no significance for predicting the low energy 
behaviour of the spectral function. 

Next, we proceed to a quantitative check of the rigid-band approximation: we set Eref = -E'0,2 
in the photoemission spectral function at half-filling, AQ,^{k, uj) and do not invert the sign of 
uj, so that we can directly compare the positions of peaks in this function and in the inverse 
photoemission spectrum for the doped ground state, A2^+{k,u!) (both spectra involve the 
single hole subspace in their final states, so that direct comparison of states is possible). 
Then, Fig. ^ confirms that in the 16 and 18-site cluster the final states for the lowest IPES 
peaks at all momenta off the (1, 1) direction (right panel) indeed belong to the single-hole 
band observed in PES at half-filling (the energies of the respective lowest peaks agree to 
10"^°t, essentially the limit of the Lanczos procedure). Fig. ^ shows that the same holds true 
in the 20-site cluster for all momenta except (0,0) and (vr,7r). In complete agreement with 
the RBA, the lowermost peaks of the IPES spectrum for the doped case are thus identified 
as the uppermost states in the PES spectrum for the undoped case. It is only at higher 
energies (~ J above the quasiparticle states) that there are states with appreciable weight in 
the IPES spectra which have vanishing or small weight in the PES spectrum at half-filling. 
The low energy physics thus should be completely consistent with rigid-band behaviour. 
It should be noted, that the above result is in strong contradiction to the 'large Fermi 
surface scenario' This would necessitate the assumption that the uppermost states of 
the completely filled, next-nearest neighbor hopping band observed in PES at half-filling 
simultaneously belong to a half-filled, nearest neighbor hopping (i.e. topologically different) 
band 'observed' for the two-hole ground state. It would moreover require the assumption 
that upon doping a full-scale transition to a topologically different band structure can occur, 
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while simultaneously the states next to the chemical potential remain unaffected and merely 
cross the Fermi level. 

Let us now turn to a discussion of the IPES final states in the 16 and 18 site cluster along the 
(1, 1) direction (left panel of Fig. H). Obviously, these states are not observed in y4o,-(fc,u;) 
so that we seem to have found 'new states', which were 'generated by doping'. The true 
explanation, however, is much simpler: for each momentum along (1, 1) the little group 
comprises the reflection by a plane along the (1,1) direction (which we denote by T^). The 
ground state at half- filling is totally symmetric, and consequently even under T^, whereas 
the two-hole ground state has d^2_y2 symmetry and hence is odd under T^. Consequently, for 
any single hole state |\Ef,^^i(A;)) the matrix element (^E'iy,i(fc)|cfc,a|\E'o,o) ((^j/,i(fc)|cfc „o-|^E'o,2)) 
is different from zero only if |\[',y_i(fc)) is even (odd) under Tm- The appearance of the 'new 
states' thus is simply the consequence of a group theoretical selection rule 0. 
We thus have to explain single-hole states with an odd parity under Tm (for the highly 
symmetric fc-points (0,0) and (vr, tt) they must have the full (ia.2_y2-symmetry). The totally 
symmetric single hole states observed in PES at half-filling can be understood in terms of 
the string picture P,p!^,|8|J5|] , where the hole is assumed to be dressed by a cloud of spin 
defects. We thus adopt the hypothesis |^ that the odd-parity states in question are similar 
in nature, but that the cloud of spin defects surrounding the hole has a nontrivial symmetry. 
Hence we define 

where j + x denotes the nearest neighbor of j in positive z-direction etc. When acting on 
the Neel state, dj^^ generates the four strings of length 1 beginning at site j. Their relative 
signs makes sure that the resulting state has (ij.2_y2-symmetry under rotations around j. A 
coherent superposition of such operators with momentum along (1, 1) consequently creates 
a state with the desired tranformation properties and we assume, that the states in question 
can be described by such a wave function. Then, for momenta along (1, 1), in Ao,_(fc, u) we 
replace Ck^a by the Fourier transform of dj^^. We denote the resulting function by y4o,_(fc, u) 
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and it is shown in Fig. ^. Again we choose Eref = £-0,2 to facihate comparison with the IPES 
spectra for the two-hole ground state. AQ^^{k, uj) is similar in character to Ao^^{k, uj) in that 
there is a 'quasiparticle peak' at the bottom of an incoherent continuum. We conclude that 
there indeed exists a band of many-body states, where a hole surrounded by a spin defect 
cloud with intrinsic (i2.2_j^2 -symmetry propagates coherently. Next, the dominant peaks in 
the IPES spectra along (1, 1) fall precisely into this band (the energies of the respective 
lowest peaks agree to an accuracy of 10~^°), so that we have clarified their nature. 
As a more quantitative check of the RBA we next consider the weight of the peaks near 
Ef] Fig. ^ compares the weight of the peak at (7r/2,7r/2) in AQ^^{k,uj) (where it equals 
the weight at (vr,0)) and A2-{k,uj) as well as the sum of the weights of the lowest peak in 
A2,+ {k,u!) and highest peak in A2-{k, —u) at (vr, 0). The RBA predicts all three quantities 
to be equal and Fig. |^ shows that they indeed agree remarkably well over a wide range of 
t/J. 

Finally, let us discuss the 'band structure' near Ef, summarized in Fig. |^ for the 16 and 
18-site cluster. In A2-{k,uj) we observe very much the same band as in AQ^_{k,uj), the 
dispersion of the peaks close to Ef being practically identical to that for the half- filled band. 
For the doped case, there is moreover an obvious correlation between the peak intensity and 
the distance from Ep as one would expect it for a Fermi liquid: comparison of Figs. || and 
1^ shows that sharp peaks exist for those momenta which are closest to the Fermi energy in 
their respective clusters, weak or diffuse peaks are seen for momenta which are more distant 
from Ep (for momenta ouside the AFBZ there is in addition a depletion of intensity over the 
whole width of the spectra, as necessitated by the kinetic-energy sum rule (||)). The same 
overall trend can also be seen in the 20-site cluster (Fig. In IPES for the 16 and 18 site 
cluster, the situation is more complicated due to the novel selection rule for momenta along 
the (1, 1) direction. Away from this line we observe in A2^+{k, uj) a portion of the band seen 
in AQ,^{k,uj) (as is the case for all momenta in the 20-site cluster, where no selection rule 
exists). In the sense of the RBA these states have partially crossed the Fermi energy; due 
to the interaction between the holes, however, there is no Fermi surface but rather a zone of 
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partially occupied momenta (indicated by the box in Fig. |^) where the quasiparticle peaks 
are split between PES and IPES. For IPES along (1, l)-direction in the 16 and 18-site cluster 
the selection rule prohibits transitions back into the single-hole band and a different band 
of many-body states with odd parity under reflection by the (1, 1) axis is seen. We have 
identified them as a hole dressed by a cloud of spin defects with intrinsic dx2_y2 symmetry. 
This band has almost no dispersion and for all momenta accessible to our diagonalization 
study remains ~2J above Ef, and thus is unrelated to any low-energy physics. 
While there is no well-defined Fermi surface in the bound states we have studied, the obvi- 
ous validity of the rigid band approximation suggests that if a Fermi surface exists at all, 
it takes the form of small hole-pockets. The precise location of these pockets in an infi- 
nite system is impossible to predict on the basis of exact diagonalization results; due to the 
near-degeneracy of the states near the surface of the magnetic Brillouin zone, the interaction 
between the holes is decisive. Whereas the 16 and 18 site cluster both suggest (vr, 0) as the 
locus of the pockets, this momentum seems to be largely unoccupied by holes in the 20-site 
cluster, where the largest shift from PES to IPES occurs at (27r/5, 47r/5). The best one 
can say is that there seems to be a trend for hole occupation at or near (vr, 0). Since it is 
the interaction between the holes, which favours these fc-points, the apparent contradiction 
with the well known fact that the single hole ground state has momentuma (7r/2, 7r/2) is not 
surprising. 

Rigid band behaviour and hole pockets in the t—J model are consistently suggested by a 
number of previous exact diagonalization works. Poilblanc and Dagotto [|1^ studied the 
A{k,Lj) for single hole states and concluded that the two-hole ground state in the 4x4 
cluster shows hole pockets at (vr,0), in agreement with the present result. Stephan and 
Horsch |I|] studied n{k) and A{k,Lj) for the two-hole ground state and concluded on the 
contrary that there is neither rigid band behaviour nor hole pockets. It should be noted, 
however, that whereas Poilblanc and Dagotto employed a quantitative criterion (presence or 
absence of a quasiparticle peak at the position of the two-hole ground state energy), Stephan 
and Horsch based their conclusions solely on the qualitative inspection of a rather limited 
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data set; as discussed above (see Tab. |) n{k) is not reliable for assigning the Fermi surface 
and our results for the spectral function in the 20-site cluster show for example that the 
Luttinger Fermi surface assigned there by Stephan and Horsch does not exist. Next, Castillo 
and Balseiro |]1T| computed the Hall constant and found its sign near half-filling to be con- 



sistent with a hole-like Fermi surface. Gooding et al. |]T^ studied the doping dependence of 
the spin correlation function in clusters with special geometry and also found indications of 
rigid-band behaviour. Finally, a study of n{k) in the range J > t (where the incoheren con- 
tinua are neghgible) with an added density repulsion to preclude hole clustering shows 
unambiguous hole pockets. It seems fair to say that the available numerical results for the 
t—J model, when interpreted with care are all consistent with rigid band behaviour and/or 
hole pockets. 

In summary, we have performed a detailed study of the doping dependence of the single 
particle spectral function up to the largest clusters that are numerically tractable. The 
results show unambiguously that rigid-band behaviour is realized in small clusters of t—J 
model: near the chemical potential, the main effect of the doping consists in moving the 
Fermi energy into the 'band' of peaks observed at half-filling. Thereby the parts of the 
quasiparticle band which remain on the photoemisson side are essentially unaffected, the 
uppermost states of this band simply cross the Fermi level and reappear as the lowermost 
states of the inverse photoemission spectrum. This behaviour is always realized, unless it is 
prohibited by a trivial reason, namely a symmetry related selection rule. In the latter case, 
there is no low energy IPES weight at all. On the PES side, modifications of the spectral 
function which deviate from the rigid band picture occur predominantly at energies far from 
Ef, and hence should be unrelated to any Fermi surface physics. In particular, the gains and 
losses of PES weight which lead to the formation of a 'large Fermi surface' in the momentum 
distribution upon doping originate from addition and depletion of incoherent weight deep 
below the Fermi energy. On the IPES side, there is no indication of the emerging of 'new 
states' at low energies in the course of doping. Such new states in the IPES spectrum are 
seen only at high energies, and hence also should be unrelated to any low energy physics. 
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Available exact diagonalization data are all in all consistent with this interpretation. 
Then, we are left with the problem to reconcile the emerging picture for the t—J model with 
experiments on high-temperature superconductors. While some transport properties are 
quite consistent with a rigid-band/hole pocket scenario ||T^, the main problem is with angu- 
lar resolved photoemission experiments ||T5[. These show peaks which disperse towards the 



Fermi energy and vanish there, as one would expect for a band crossing. Thereby the locus 
of the 'crossing points' is remarkably consistent with the predictions of band theory, which 
in turn is inconsistent with hole pockets. In a Fermi liquid, the contradicting quantities 
actually fall into two distinct classes: photoemission spectra depend on the wave function 
renormalization constant Z^, transport properties do not. Hence, if one wants to resolve the 
discrepancy concerning the volume of the Fermi surface entirely within a Fermi liquid-like 
picture, the simplest way would be to assume a 'small' Fermi surface (to model the trans- 
port properties) and explain the photoemission results by a systematic variation of Zh along 



the band which forms the Fermi surface, i.e. similar to the 'shadow band' picture [|T^. A 
trivial argument for such a strong variation in Z^ would be that irrespective of the actual 
band structure, a distribution of PES weight in the Brillouin zone that resembles the noin- 
teracting (band theory) Fermi surface always optimizes the expectation value of the kinetic 
energy. A wave function which gives a substantially different distribution of spectral weight 
consequently has a very unfavourable kinetic energy and hence is ruled as the ground state 
from the very beginning. Also, it should be noted that the spectral weight in the shadow 
band should even decrease when the charge fluctuations in the original Hubbard model are 
taken into account: when going back to the Hubbard model, the (negative) exchange energy 
in the t—J model then is split into a (positive) contribution from the Hubbard repulsion and 
a (negative) gain in kinetic energy, so that the total expectation value of the kinetic energy 
certainly becomes more negative than for the t—J model. Via the kinetic-energy sum rule 
(0), one can infer that this necessarily leads to an even more free-electron like distribution 
of spectral weight. All in all, adopting the rigid-band behaviour found above, the available 
diagonalization data for the t—J model [|li[To|-[T^ as a whole then are reasonably consistent 
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with the above scenario. 
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TABLES 

TABLE I. Momentum distribution in the ground state of of the t—J^ model {t/Jz = 2) in 
a staggered magnetic field (O.lJz) and the t—J model (t/J=2) in a 4x4 cluster. To stabilize 
a ground state with the same point group symmetry (Bi) as for the t—J model, a 2"*^ nearest 

ncigiilior h()p]")iug lorm of sircuglli —. 7-^/10 has V^ccu added lo llic I — J- UK)d(^l. 



k 


(0,0) 


(f>0) 


(7r,0) 


V2' 2/ 




(7r,7r) 


n{k){t- 


Jz) 


0.5481 


0.5374 


0.3491 


0.4993 


0.3387 


0.2522 


n{k){t- 


J) 


0.5565 


0.5421 


0.3378 


0.4974 


0.3225 


0.3197 
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FIGURES 

FIG. 1. Comparison of the photoemsission spectra at half- filling and in the ground state with 
2 holes: AQ^-{k, —u) (dotted line) A2-{k,—Lj) (full line) A2^+{k,uj) (dashed-dotted line) for the 
16 and 18-site cluster. 5 functions have been replaced by Lorentzians of width O.lt. 

FIG. 2. Same as Fig. |I| but for the 20-site cluster. 

FIG. 3. PES spectrum Ao^_(fe,Lj) (full line) and IPES spectrum A2^+{k,uj) (dashed-dotted 
line) for the 16 and 18-site cluster. The reference energies are identical so that direct comparison 
of the peaks is possible, 6 functions have been replaced by Lorentzians of width O.OSt. 

FIG. 4. Same as Fig. ^ but for the 20-site cluster. 

FIG. 5. PES spectrum AQ^^{k,uj) (full line) and IPES spectrum A2^+{k,uj) (dashed-dotted 
line) for momenta along (1, 1) in the 16 and 18-site cluster. The reference energies are identical so 
that direct comparison of the peaks is possible, 5 functions have been replaced by Lorentzians of 
width 0.05t. 

FIG. 6. Comparison of the t/ J-dependence of the PES pole strength at (7r/2, 7r/2) at half-filling 
(squares) and in the two-hole ground state (triangles) and the added weights of the lowest IPES 
peak and highest PES peak at (vr, 0) in the two-hole ground state (circles) 

FIG. 7. Schematic quasiparticle band structure in the neighborhood of Ef for the 16 and 18-site 
cluster. Up-triangles (squares) give the position of the highest peak in A2^-{k, —uj) down triangles 
(circles) the position of the lowest peak in A2^+{k,uj) for the 4x4 (18-site) cluster. The positions 
of the highest peaks in AQ-{k, —lo) (dots) are also given. The various groups of spectra have been 
shifted so that the energies of the respective PES peak at (0, 0) coincide (shift between doped 16 
and 18-site cluster: 0.275i). 
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